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Abstract 

The excitation of electromagnetic (EM) wake waves in electron plasma by an 
one-dimensional bunch of charged particles has been considered in the presence of 
intense monochromatic circularly polarized electromagnetic (CPEM) pump wave. 
In the zero state (in the absence of bunch) the interaction of the pump wave 
with plasma is described by means of Maxwell equations and relativistic nonlin- 
ear hydrodynamic equations of cold plasma. The excitation of linear waves by 
one-dimensional bunch is considered on this background. It is shown that there 
are three types of solutions of linear equations obtained for induced waves corre- 
sponding to three ranges of parameter values of the pump wave, bunch and plasma. 
In the first range of parameter values the amplitude of transverse components of 
induced waves is shown to grow as the bunch energy and after some value of the 
relativistic factor of the bunch to be almost independent of the energy and increase 
proportional to the intensity and frequency of the pump wave. The dependence 
of longitudinal component of induced waves on the relativistic factor of the bunch 
is weak. Its amplitude and wavelength grow as the intensity of pump wave. The 
second range of parameters is a resonance one. The amplitude of the wave excited 
by the bunch is a linear function of the distance to the bunch. In the third range 
of parameter values the longitudinal component of induced fields are localized near 
the bunch boundaries and are exponentially decreased with the increase in distance 
from these boundaries. The amplitudes of transverse components of induced waves 
reach a constant value with the distance from the bunch boundaries. 
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1 Introduction 



At present the studies on new methods for charged particle acceleration by means of 
wake fields generated in plasma by laser radiation (BWA (Beat Wave Acceleration), 
LWFA (Laser Wake Field Acceleration)) and by bunches of relativistic particles 
(PWFA (Plasma Wake-Field Acceleration)) in flight through plasma are intensively 
developed (see, e.g., the reviews [f, 2] and references therein). The intensity of 
acceleration fields (in the order of 10 7 — W 8 V/ cm), attained by these methods can 
be used both for the charge acceleration, and for focusing of electron (positron) 
bunches in order to obtain the beams of high density and to ensure high luminosity 
in linear colliders of next generation. 

The linear theory of wake field generation by rigid bunches of charged particles 
in boundless and limited plasma was developed in many works (see, e.g., [3-11]). 
The nonlinear theory of wake field generation by a rigid one-dimensional bunch 
of finite extent was developed in [12-18]. An important result of this theory is 
the proof that the wave breaking limit is E max = (mv Lu p /e) [2(-f — l)] 1 ^ 2 and 

it is reached at rib/no < 1/(2 + Jq 1 ), where lu p = {AimQe 2 /m) 1/>2 is the plasma 
frequency of electrons, n& and no are the densities of the bunch and plasma electrons, 

1/2 

7o = (l — (3q) is the relativistic factor of the bunch, /3 = v /c. The Dawson 
[17] wave breaking limit is equal to E mayi w 2mvQL0 p /e when 70 ~ 1 (/?o <C 1). In 
the linear case n^/no <C<g; 1, E max ~ 2(mvQbj. p /e)(rn,/nv) for arbitrary 70. 

The one-dimensional relativistic strong waves can be excited in the plasma by 
wide relativistic bunches of charged particles or intense laser pulses [1, 19], (when 
k p ao ^> 1, where k p = lo p /vq, ciq are the characteristic transverse sizes of bunches 
or pulses). 

In the present work the effect of the CPEM pump wave with an arbitrary in- 
tensity parameter (A = eE^/mcuj^, where E and too are the amplitude and the 
frequency of the EM wave) on the excitation of EM wake waves by one-dimensional 
relativistic electron bunch in cold plasma has been studied. Here the rate of elec- 
tron oscillations in the pump wave may be of the order of light velocity. One can 
obtain the exact solutions of the Maxwell equations and of nonlinear hydrodynamic 
equations [18, 20] for CPEM wave interacting with plasma as well as derive the ex- 
act dispersion equation for waves propagating in the same direction as pump wave 
[20, 21]. The parametric instability of the plasma in the presence of CPEM wave 
has been studied rather well yet in early works (see, e.g., [20, 22, 23] and literature 
therein). Max and Perkins [22] investigated an aperiodic low frequency instability 
of plasma in the dipole approximation. The instability of plasma in the presence 
of strong CPEM wave was considered in [20] and it was shown that at parametric 
excitation of nonpotcntial oscillations in plasma by a CPEM the relativistic motion 
of electrons is essential for the arbitrary value of pump wave amplitude. 

This paper is organized as follows. The derivation of equations for the velocity 
of motion of an electron liquid in plasma and for EM fields excited by an one- 
dimensional electron bunch in the presence of CPEM wave is given in Sec. II. The 
interaction of the pump wave with plasma (in the absence of a bunch) is described 
by Maxwell equations and nonlinear hydrodynamic equations of cold plasma. In 
this case a spatially homogeneous state of plasma is possible [18, 20]. Then, assum- 
ing that this state is weakly perturbed by one-dimensional bunch in the plasma, the 
equations for induced fields, density and velocity of plasma electrons are obtained 
using the perturbation theory methods. In Sec. Ill the general expressions for 
induced EM fields of one-dimensional bunch with an arbitrary density profile are 
obtained using the method of Green's functions. There are three ranges of values 
of plasma, pump wave and bunch parameters, where the properties of the Green's 
functions and, hence, the behavior of excited waves, are abruptly changed. In Sec. 
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IV the fields excited in plasma by a bunch with uniform distribution of electron 
density are considered. For each of the mentioned ranges of parameters the ex- 
pressions for induced fields components are obtained and corresponding numerical 
calculations are carried out. A brief summary of the results is given in Sec. V. An 
alternative method for obtaining the induced fields is discussed in the Appendix. 



2 Basic Equations 



As an initial system of equations we shall use the Maxwell equations and relativistic 
hydrodynamic equations of motion of cold electron plasma under assumption that 
the oscillation velocity of plasma electrons in a CPEM wave much exceeds their 
thermal velocities and the frequency coo of the CPEM wave is much higher than the 
frequency of electron-ion collisions: 



V x B = 



1 <9E 

c~dt 



Aire Aire liT 
nv uN b £ , 

c c 



V x E = — 



1 dB 



VB = 0, 



(1) 



(2) 



VE = -47re (n - n ) - 47reiY b (£), 



(3) 



— + vV v = 

ot i 



1 v 
E + -v x B - — (vE) 



- + V(nv) 



0, 



(4) 



(5) 



where no is the unperturbed electron density, Nb(£) is the density of one-dimensional 
electron bunch moving with velocity u = ue z , £ = z — ut. 

In the CPEM wave propagating along the z axis a spatially homogeneous state 
of plasma may be established, in which the EM field and velocity of electrons will 
be determined by expressions [18, 20] 



E = E (e x cos C + e y sin () , E { 



0z 



0, 



(6) 



B = —E (-e x sin ( + e y cos Q , 



B 0z = 0, 



(7) 



v e = cf3 e (-e x sin ( + e y cos () 



0. 



(8) 



where C = u t - k z, k = {lj /c)^s{lo ), e{uj) = 1 - ujI/lj 2 , lo\ = - A 

P e =v e /c, A = eEo/mcLJo, 



A 



v e = c 



\f\TA 1 ' 



(9) 



4nnoe 2 /m is the plasma frequency, c is the velocity of light. 



Consider small perturbations of plasma due to the presence of an electron bunch 
with density of Nb. The linearization implies that the analysis is valid only on 
condition that the bunch density Nb -C no. We shall write all quantities in the form 
/ = fo + f , where /o are determined by means of Eqs. (6)-(9). The linearization of 
the system of Eqs. (l)-(5) gives the following set of equations for induced variables 
/' 
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V x B = 



1 <9E 47re 



c dt 



„ - ldB 

VxE =—- s", 

c at 



VE 



ft i \ 47re 
^n v + n v e J — u^Vfe(^), 

VB' = 0, 

4iren - 4ireN b (£), 



(10) 

(11) 
(12) 



dv 

JC_ + (v e V)v' + (v'v) v e 



+ - 



("•) 



E' + -v e xB' + VxBo-^ (v e E' 



E + -v e x B 



v E f 



(13) 



+ 



<9n 

~dt 



n Vv + (v e V) n = 0. 



(14) 



The set of Eqs. (10)-(14) is a system of partial differential equations with periodical 
coefficients with respect to the variable (. 

If we pass from x and y components of fields and velocities of plasma electrons 
to new variables 




E x ± iE y 
B' x ± W y 




(15) 



this will mean that we convert these into the rotating system of reference of the 
CPEM wave. Here the set of Eqs. (10)-(14) passes into a system of ordinary 
inhomogencous differential equations with constant coefficients 



1 8 2 £ ± 

7 2 de 

4iren 



T 2i (ko - 0- 



,2 _ ^0 



dV 



-r- -^vi , a n i0 e dn 
<9£ c n n d£ 



(16) 



1 8 2 B ± 



7 2 ^2 



2i (k - 



.2 ^0 

fc ° C 2 



6 d 



± 



, , . , n . iv e dn 

T ikoV* + k v e — ± — 

<9? 04 



(17) 



2v± 



a 2 v 



0d± 



d£ 2 



— Vwf 

1 — — 



a£ u ) 2 



(18) 



-Aireri - 4neN b (£), B' z = 0, 



(19) 



4 



a£ mu 



kg Eg 

2uj 



(v + + v-) 



(20) 



no — , 

u 



(21) 



where f3 = u/c. At the derivation of Eqs. (16)-(21) we assumed that all quantities 
depended on the variable £ = z — ut. One can see from these equations that owing 
to the presence of CPEM wave in the plasma, the induced fields and velocities of 
electron liquid motion are parametrically related. Besides the longitudinal compo- 
nent of the induced electric field the right-hand side of Eq. (20) for the induced 
velocity longitudinal component comprises also the Lorentz force. The second term 
arises on account of the interaction of unperturbed velocity of plasma electrons with 
the induced magnetic field. The third term is due to the interaction of induced ve- 
locities of electrons with the magnetic field of CPEM wave. As is seen from the 
expression (20), the coefficient of the second term (j3 e ) in the square brackets tends 
to 1/2 and the coefficient of the last term grows with the intensity of the CPEM 
wave. Below we shall examine the solutions of Eqs. (16)-(21). 



3 Calculations of Green's Functions 

To solve the system of Eqs. (16)-(21) one can obtain the equations for each of 
variables n , E z , V , B ± 1 £ ± . Such an approach to determination of induced 
fields and velocities of electrons has its advantages and is discussed in short in the 
Appendix. Here we shall solve the system of Eqs. (16)-(21) by expansion of induced 
quantities in Fourier integral in variable £. Then, after some transformations we find 

/ Gi e) (£'-0 \ 
dS'N b {Z')\ G^ e) (£'-£) , (22) 

where G { z e) ,G { l > , G ( " l) are the Green's functions respectively for quantities E z , £ + 
and B + that are as follows: 




/+oo 
-OO 



dkDi(k,uj) iks 



k D(k,u) 

+ °° 



Gy(s) = / dk(ku + Lu Q ) 

c J — oo 

/+oo 
dk(k + k 



D{k,u>) 



MM it. 

D(k,uj) 



where u = ku, 7 2 ~ 1 — [3 2 . Here we made the following notations: 



(23) 
(24) 
(25) 



R±i(k,u>) = (k±k ) 2 - ^-^-^-e{u) ± u ), 



(26) 



a2 2 

D!(k,uj) = uj 2 R 1 (k,uj)R- 1 (k,u;) + [ k 2 - 



u> 



[R 1 {k,u,) + R- 1 {k,u ) )}, (27) 



2, ,2 r 



a2 

D(k,w) = Lu 2 e(uj)R 1 (k,uj)R^ 1 (k,uj) + ^— ^ 



[R 1 (k,Lu) + R- 1 (k,u>)]. 

(28) 



5 



From expressions (15) one can find the transverse components of induced electric 
and magnetic fields, as well as the velocities of plasma electron motion by taking 
the real or imaginary parts from complex quantities £+, B + and w + . As a result 
we have 

jfe! w *$wr saw <»> 

where 



B'(z,t) J~ \ mo ;-^-v ^(0 



{ Br® )- Rc { B+® ) ' { B*® )- lm { B+® ) ■ (31) 

Thus, one can see from expressions (29) and (30) that the transverse compo- 
nents of induced fields describe the modulated oscillations in plasma. E.g., one can 
represent E y in the form 

tf y (z,t) = E ±Q ®sm(C + M€))- (32) 

Here E± ® = y/E?®+Ef® is the amplitude and Vo(£) = axctan[Ei® / E r ®] 
is the phase displacement of oscillations (other induced quantities may be similarly 
written in the form of (32)). In the system of reference of bunch rest the expression 
(32) describes the transverse harmonic wave. In an arbitrary reference system the 
expression (32) describes a modulated transverse wave, the profile of which is given 
by the function E±_o(Q- Note also that in the absence of CPEM wave (J3 e = 0) 
all transverse components of induced quantities turn to zero, and Eq. (19) and its 
solution (the first expression in Eq. (22) together with expressions (27) and (28)) 
turn into the known expressions for linear one-dimensional fields. 

Now calculate the Green's functions that are determined by expressions (23)- 
(25). The poles of integrals in the expressions (23)-(25) are determined by the 

dispersion equation D(k,uj) = 0, that has been comprehensively studied for an 

v 

arbitrary case in Ref. [20] (i.e., without the Cherenkov condition lo = ku to be im- 
posed). In the absence of a CPEM wave (/3 e = 0) we obtain dispersion equations for 
ordinary plasma and transverse (EM) waves from the expression (28) with to = lu p 
and uj 2 = u)p + k 2 c 2 respectively. In the presence of a pump wave ((3 e ^ 0) there 
arise coupling waves in the plasma, the growth increment of which for small ampli- 
tudes of the CPEM wave (f3 e <C 1) linearly increases as j3 e . Thus, the parametric 
instability of coupling waves persists down to the value of (3 e = 0. However, in the 
case of non-dense plasma (loq < 10 13 sec -1 or no < 10 17 cm~ 3 ) and for optical values 
of loq (u>q ~ 10 15 sec -1 ), the time of relativistic bunch interaction is much less than 
the time of parametric instability development [20] and here we shall not make any 

allowance for their effect on the excitation of waves in plasma. 

v 

In case when the Cherenkov condition is approached (uj = ku) we obtain from 
expression (28) the following dispersion equation 




If we introduce a dimensionless wave vector A that is related to k by means of 
expression k — (u)l/u)\, then the solution of Eq. (3) will take on the form 
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where 

F = v^A 2 - 1 - /3aA, (35) 

A = loo/uj p , a 2 = Vl + A 2 . 

The behavior of solutions of Eq. (33) and, consequently, the nature of induced 
fields and of velocity of plasma electrons strongly depend on the sign of radicand 
in the Eq. (34). In the following we shall identify the regions where the radicand 
is positive, equal to zero or negative as the I st , II nd , and III rd regions respectively. 
First, we shall determine the boundaries of the region where this radicand turns 
to zero. If the radicand is made equal to zero one finds for the boundary of the 
above-mentioned region 



A=^( Jl + 1 2 F 2 ±p 1 F + ), (36) 



when a > 1, and 



A=^( A /l + 7^ ±/ 3 7 F_), (37) 



when 1 < a < ao(j), where 



. 1/4 



The Eqs. (36)-(39) were obtained on the assumption that 7 > 71 = 1.45, where 71 
is a real positive root of equation 27 2 (7 2 — 2) = 7 — 1 that satisfies the condition 
7 > 1. 

For small pump wave amplitude (a — 1 <C 1) the values of functions given by 
the first expressions (with plus sign) and those of functions given by the second 
expressions (with minus sign) of Eqs. (36) and (37) coincide. We denote these 
values respectively as A± and from Eqs. (36) and (37) we can obtain for them the 
following expressions: 



A± =Y + 4^=i) ± 5- (40) 

In Figs. 1,2, and 3 the regions I, II and III (the lines coincide with the region 
II) are given for values of 7 = 1.5, 7 = 20 and 7 = 100 respectively. The lines in 
these figures are closed at the infinity (i.e., when a — > 00). For modern lasers with 
intensities II < 10 20 W/cm 2 we have the following restriction on the parameter a: 
a < 3 for values of frequency of about uj — 3 x 10 15 sec _1 . So, as is seen from 
Figs. 1, 2 and 3 and Eqs. (36) and (37), for values of parameters luq ~ 10 15 sec _1 , 
n < 10 17 cm~ 3 (oj p < 10 13 sec~ 1 ) the solutions of Eq. (35) are in the region I for 
wide spread of the values of 7 (up to the values of 7 < (uj /oj p )a and higher). The 
solutions of Eq. (33) will be in the regions III or II if the condition 7 ~ aA = 
(bj /Lo p )a > 100 is observed (the ultrarelativistic bunch). 

Now we shall calculate the Green's function for each of the regions in separate. 
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3.1 Region I 



In this case the roots of Eq. (33) are real and lie within the upper half plane of 
the complex variable k. Integrating the expressions (23)- (25) over the variable k we 
find 



4ire6(s) 
\\ - X 2 _ 

X 2 _ + 



4/3 2 7 4 F 2 ) cos ( '--k p s 



(41) 



me 



AF 



F 2 



4a 4 ' 



[0(a)-B(-a)] + 



. . , 2/3 7 2 aAF\ /A_ , 

+1 ( 1 ^2 I C0S I S S 



1 



■ sin 



A + \ (/ 

,2 



A.. 



2/3 7 2 aAF\ 



1 






sin [ — k 


a7 


V a 


*± 








a 



G ( , m) (s) 



= — 7re- 



V^ 4 



FFn 



F 2 + ^l- 2 



W) 0{-s)\ + ^S-V ) x (43) 



2[3 2 -/ 2 FF 



-0 (2 7 2 F - F ) 



cos 



: fc p s 



- 1 



A 2 . - A 2 

2/3 2 7 2 F J Fi,\ / A 

^2 I C0S 



■ kpS 



-L sinf^vVf sinf^V 
A + \ a / A_ \ a 



where fe p = a> p /u, Fo = \/a 2 A 2 — 1, 0(s) is the Heavyside function. It is seen from 
Eqs. (41)-(43) that the bunch excites in plasma the oscillations of two types with 
frequencies ujlX+ and ojlX-. In the a = 1 limit (the CPEM wave is absent) 
and G^ turn to zero and for G^ we obtain the expression 



Gi e) (s) = 47re(9(s)cos(fcps), 
that coincides with those given in Refs. [24, 25]. 



(44) 



3.2 Region II 

In the second region the radicand in Eq. (34) turns zero. In this case A 2 ^ = X 2 _ = A 2 ,, 



A 2 = ^ + 2/? 2 7 4 F 2 = 1 ± 7 



1 



>o, 



(45) 



where the sign "+" corresponds to the values of parameters 1 < a, 7 > 7 i (see the 
Eq. (36)). The sign "-" corresponds to the values of parameters 1 < a < ao( 7 ), 
7 > 7 i (see the Eq. (37)). The poles ±Ao in this region are multiple ones and lie in 
the upper complex plane k. The calculation of integrals in Eqs. (23)-(25) gives 



G^(s) = Airee(s) 



cos 



2a Xq 



Ao 



kpS 



(46) 



F 2 _i_ a 4 -! ., -2 



27re/3 7 2 Va 4 - 1 



6{ S ) 



A 3 , a 

(oA - 2/3 7 2 .F) — (fc p s) ^- 

a Aq 



«A - 23j 2 F + i (2/?7 2 AF 
Ac 



\ 2 
A 



cos i —kpS 



(k p s) 



sin 



-7re- 



F-Fn 



F 2 



ai-l -2 



[«(») - «(-.)] + 7 V <(«X«8) 



Aq a 



a 



+ 



^(2^ 2 7 2 FFo-A 2 )+z/3(2 7 2 F-F ) 



sin ( —/.-...-• 



Note, that one could obtain the Eqs. (46)-(48) from expressions (41)-(43) by tending 
A + — > A_ — > Ao in the latter. The uncertainty that arises here may be eliminated 
by using the L'Hopital rule. 

Unlike the Eqs. (41)-(43) the Eqs. (46)-(48) contain oscillating terms, the 
amplitudes of which are linear functions of s. We shall see below that such a 
dependence of Green's functions leads to an exponential growth of the induced 
fields as a function of coordinates. 

3.3 Region III 

In this case the roots of dispersion Eq. (33) are complex quantities and are dis- 
tributed over the complex plane A symmetrical with respect to the origin of coor- 
dinates. We shall denote the roots in the upper half plane as A± = ±a + i5 (the 
roots in the lower complex half plane will be ±a — iS), where 



a 




^JP 2 + P 2 ± P l 



P 1 = ^ + 2/?VF 2 , 



(49) 
(50) 



//3 2 7 2 



a 4 - 1 



4 7 2 F 2 _ p2 



(51) 



After calculation of the integrals over k in Eqs. (23)- (25) we have the following 
expressions for the Green's functions: 



ne [6(s) — 0(—s)] exp 



■-fep|s| 



1 

aS 



— la 2 -5 2 + 



1 



4f3 2 j 4 F 2 I sin 



(52) 

+ 2 cos (^—kpS^j 



G ( l\s) = me 



AF 



F 2 



[0(8)-6{-8) 



1 



aS (a 2 + 5 2 ) a 2 



■ exp 
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x < a 



(oA - 2/3^F) i±l + 



Ai8d-i 2 aAF 



COS ( — fcr,S 

a 



+ 



sin 



GV m) (s) 



Va 4 - 1 FF 



a2 F2 + ^7- 2 



[*(«)-*(-*)] 



■exp fep|s| J5A) 



a 6 (a 2 + S 2 ) a 2 



4<5/3 7 2 FF \s\ 



+ 



a 2 + S 2 

\ 2 



i (2 7 2 F - F ) 



cos 



+ 



(a 2 + 5 2 y - 2(3 2 1 2 FF (a 2 - S 2 ) \ s \ 
a 2 + 5 2 s 



+ 05 (2j 2 F - F ) 



sin(^ p | S |) 



It follows from Eqs. (52)- (54) that in the region III the Green's functions contain 
terms that oscillate with exponentially decreasing amplitude. 



4 Induced Fields for Specific Choice of iV&(£) 

In this section we shall calculate and investigate the induced fields for a specific 
profile of electron bunch density. We shall assume that the electrons are homoge- 
neously distributed in the bunch with the density n b (n&<no), and the length of 
bunch is d, i.e., 



N b (0 = n b [9(0 - 6(£ - d)} . (55) 

As was noted in Sec. Ill, the values of induced fields strongly depend on the 
fact, in which of three mentioned regions lie the values of the discriminant of Eq. 
(33). The induced fields in each region (I, II or III) of the values of bunch, plasma 
and pump wave parameters will be considered separately. 



4.1 Region I 

In this case the Green's functions are determined by the Eqs. (41)-(43). The 
substitution of these expressions and Eq. (55) into the Eq. (22) will result in the 
obtaining of following expressions for the induced fields ahead of (£ > d), inside 
(0 < £ < d) and behind (£ < 0) the bunch: 

£ >d; 

<(0=0, (56) 



£+(0 = -iE — {k p d) 
no 



Va 4- ! 



BAF 



n P ' a 2 4F 2 + 4f i 7" 



4F 2 + ^ T V 2: 
pFF 



(57) 
(58) 



< £ < d; 



-^( A - + ^- 4 ^ 2 ) sta (¥^-«)}- 
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£+(0 



iE, 



n b /3\/a 4 - 1 



no a 



AF 
4F 2 + 



fcp(d-20- 



07 



2 r \2 



A 2 , - A 2 . 



A 2 -2/3 7 2 aAF . /A+ 



4 



sin I — k p (d — £) 



(60) 



\ 2 _-2f3 1 2 aAF . /A 

— sin I —k p (d - £) J + 



+i (aA - 2/? 7 2 F) 



^(l-cos(^Md-O))- 



A 3 



U - cos f^-fc p (d-£) 



C<0; 



n 6 



/3Va 4 - 1 



FFn 



no or 



~\\-\ 2 _ 



4F 2 + ^-j- 2 



k p {d-2£)+ 



(61) 



rA^- 2/3 2 7 2 FF . /A+ 



sin - 



XL - 2/3VF F . fx 
A 3 _ 



sin — k p (d-£) - 



i/3 (2 1 2 F - F ) 



^ (i - - - fl)) - -r (i - - - 0))] ] } • 



^(0 = 



ni, af3 
no X\ — A 2 . [ A^ 



+ 

A 2 _ 



xl + a - 



- 4/3 2 7 4 F 2 



4 - 1 



sin ( —k p (d — £) 



sin^Ap(d-o) +sin^±Ap^|63) 

A_ 



£ + (0 = 



/Va 4 - 1 



AF 



ng a 



4F 2 + ^V 2 
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A 2 , - A 2 



A 2 - 2(3j 2 aAF 
4 



(fcpd)- 



(63) 



sin I — -k p (d ~ £)) + s i n ( ~^p? 
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A 3 

+i (aA - 2/?7 2 F) 
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-jj- ( cos ( — kpS, ) - cos ( —k p (d - £) 
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+ 
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\\ - A 2 
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Ai 

(2 1 2 F - Fo) 



sin [ k p (d — £) ] + sin [ — fc p £ 



A^~ (° OS ("a^^) ~ C ° S (~<t kp( " d ~ ^ 



A 2 



cos 



A_ 



k p £ - cos k p (d - £) 



Here we have adopted the notation Eq = mcui p /e. 
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4.2 Region II 

In this range of values of bunch, plasma and pump wave parameters the induced 
fields ahead of the bunch (£ > d) coincide with the values of induced EM fields in 
the region I and are given by the expressions (56)- (58). The induced fields inside 
and behind the bunch are 

< £ < d; 



n ZA 



A \ a J \ a 

(65) 



iE 



n b [3^a A - 1 



AF 



n 
2Ag 



k p (d - 20 + ^ (aA - 2/3 7 2 ^) - (66) 
o 



+A 



6/3 7 2 aAF - Aq + iX 2 (aA - 2/3 7 2 F) 
(Aq - 2/3 7 2 aAF) ^fcll + 2i (aA - 2/3 7 2 F) 



sin ( —k p (d - £) ] -+ 



cos I —kJd — £) 
1 a 



Ti 6 /Va 4 - 1 
n a 2 



M<* - 2 + ^5 [ 2 ^ A o ( 2 T 2j7 - ^o) 4(67) 



4p2 _|_ a 4 -U -- 2 



+ 6(3 2 j 2 FF - A 2 - z/3A 2 (2 7 2 F - F ) 
+A 



k p (d-0 



sin | —k p (d - £) ) + 



(A 2 - 2/5 2 7 2 ^F ) fcp(d - 2z/3 (2 7 2 F - F ) 



< os I — fc p (d — 



£<0; 



^(0 = ^ 



n 2A 2 , 
+(1 - Ag) 



3A 2 



Ao 



sin 



k p{d-Cl ) +sinl — 



//,,('/ .';) COS ( ^fc p (d - £) j + {kpC) COS f "^pC 



(68) 



£ + (0 



iE 



n b (3Va 4 - 1 [ AFk p d 



n a 

ML 

2X1 



(69) 



6/3 7 2 aAF — Aq + iX 2 (aA - 2/3 7 2 ^ 



k P (d-0 



Ao 



sin | —k p (d - £) ) + 



+ 
+A 
+A 



6/3 7 2 aAF - A 2 - U 2 (aA - 2/? 7 2 F) M 



sin ( — k p £, 
1 a 



(Aq - 2/3 7 2 aAF) + 2 i (aA - 2/3 7 2 ^) 



cos — k p (d — £)) + 



(A 2 - 2/3 7 2 aAF) M - 2* (aA - 2/3 7 2 F) 



o is ( —/,',,£ 



-Bo 



rib /3\/a 



FF 



n a z 



4F 2 + a 1 -- 2 



{k P d)+ 



(70) 
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a 7 
2Ag 



+ 
+A 
+A 



6/3 2 7 2 FF - A 2 , — z/3A 2 (2 7 2 F - F„) 
6/3 2 7 2 FF — Ag + z/3A 2 (2 7 2 F - F ) M 



sin ( ^-k p (d - f) 



sin ( ^A>£ ) r 



(A 2 -2/? 2 7 2 FF )Mr__J) 



2ip (2 7 2 F - F ) 



cos ( ^-k p (d - £) ) + 



(A 2 - 2/3 2 7 2 FF ) M + 2t/3 (2 7 2 F - F ) 



COS ( — 



where Aq is determined by the Eq. (45). 



4.3 Region III 



In this region the Green's functions are determined by Eqs. (52)-(54). The expres- 
sion for longitudinal electric fields over the space (— oo < £ < +oo) is obtained after 
the substitution of Eqs. (52)-(54) into the Eq. (22) and appropriate integrations: 



EM = E Q 



rib a/3 
~oA{a 2 + S 2 f 
S 



--[a 2 + S 2 - 
a 



a 4 -l 



+ 4/3 2 7 4 F 2 ] x 



(71) 



+ 



+ i (a 2 + 5 2 + - 4/3 2 7 4 F 2 ] x 



exp ( — fc p |£| j cos (—kp^j — exp ( k p \£ — d\j cos ^— fc p (£ — d) 



The expressions for induced transverse electromagnetic fields ahead of (£ > d), 
inside (0 ^ £ ^ d) and behind (£ < 0) the bunch are obtained in similar manner: 

£>d; 



*+(0 - -^o^(M)^^ ^ 



n 6 



no4a(5(a 2 + (5 2 ) / a 
2c«5 (aA - 2/3 1 2 F) - (V + S 2 ) 2 + 2/3 7 2 FaA (3<5 2 - a 2 )) 

-p(-£we-*)-(=l«- fl )-^(-M-(=M) 



x (72) 



+ 



+ 



(aA - 2/3 7 2 F) (a 2 - 5 2 ) + -^-^ ((a 2 + 5 2 f - 2/3 7 2 FaA (3a 2 - 5 2 )) 



B+(0 = 



£o — (fcpd) ^ — 

n a" 1 



f3FF Q 



AF 2 



2°__l -2 



V^ 43 ! /3 7 2 



a 4a<5(a 2 + 5 2 )' 



x(73) 



x < -a 



(a 2 + <5 2 ) + 2/3 2 7 2 (3<5 2 - a 2 ) FF Q 
a 2 +S 2 



+ 2if3S (F - 2 1 2 F) 



exp ^— — k p (£ — d)^j cos ^— fc p (£ — d)^ — exp ^— -fc p £^ cos (—k p (^ 



+ 
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+ 



s£ <e s= rf; 



(a 2 + S 2 ) 2 -2f3 2 1 2 (3a 2 -5 2 )FF 2W 
S ~ " ^2 — + if* ( S " « 2 ) ( F o - 2 7 2 ^) 

cxp (—-k p t;j sin ^— k p t^J — cxp (— —k p (^ — d)j sin ^ — fc p (£ — d) 



K J °n KpK qj} a 4^2 + 4^7-2 n 4aS (a 2 + 5 2 f 



~ n b /3 2 7 2 V^T 
£jq - x 



(74) 



^2a6 (aA - 2(3j 2 F) cxp ^--fc p ^ cos (-fcp^) - cxp (--k p (d- £)j cos (-fc p (d - f) 

cxp ^— — fc p ^ sin (—k p l^J — cxp ^— — k p (d — sin ^—k p (d — g)j 



+ (aA - 2/3 7 2 f ) (<5 2 - a 2 ) 



* 2 + <S 2 



( a 2 + < 5 2 ) 2 -2/3 7 2 J FaA (a 2 - 35 2 ) 



exp ^ & P £J cos ^— fcp^ + cxp y k p (d — £)J cos ^— fc p (d — — 2 

r (a 2 +«5 2 ) 2 -2/3 7 2 faA(3 a 2 - ( 5 2 ) 



i 2 + S 2 



cxp sin (~ k pt) + CX P - sin ~ ^) } ' 



~ n b ^FFq ~ 

-i?o 5 (kJd - 2£)) ^- h S 

n a 2 V pV 4F 2 + ai_i 7 -2 % a 4aS (a 2 + S 2 Y 



x < -a 



(a 2 + S 2 ) 2 + 2/3 2 7 2 (3<5 2 - a 2 ) FF Q 
1 £ff 2 } —± + M/W (F 2 7 2 F) 





a 2 + S 2 




{a 2 + S 2 ) 2 


- 2/3 2 7 2 (3a 2 


- S 2 ) FF 


— 


a 2 + 5 2 




,(- 2 + ^ 2 ) 2 


- 2/3 2 7 2 (3a 2 


- S 2 ) FF 



(a 2 + 5 2 ) 2 + 2(3 2 j 2 (3S 2 -a 2 )FF a 
± '- — — i '- - 2if3S (F - 2j 2 F 

■+0 (S 2 -a 2 ) (F -2 1 2 F 

- i/3 (5 2 - a 2 ) (F - 2 1 2 F 



— cxp ^ k p £j cos (—kpt^J 

1 — cxp ^— -k p (d — cos (^—k p (d — £) 



£ <0; 

5+(0 = 



a 2 + S 2 



no 



a 
<5 



cxp ( — k p (d — £) j sin ^ — fc p (d - 



n b (3 2 j 2 Va^l 



x < a 



a 4F 2 + ^~f- 2 " n AaS (a 2 + S 2 ) 2 a 

2 



(*>) 



+ 



j o \ (a 2 +£ 2 ) 2 -2/? 7 2 FaA(a 2 -3<5 2 ) 
25 (aA - 2/3 7 2 F) + 1 J* +s2 \ L 

;,(£«)„(=«).,(-£ Mi - <) ).(= wl .o 

2 \ / 2 c-2\ r 2/37 2 J FaA(3a 2 -(5 2 )-(a 2 + ( 5 2 ) 2 
(aA - 2/?7 2 f) (^ - S 2 ) + l 8^1 1 ~ 2 + / 2 [ >- 

cxp (J^kp^j sin {j^k p ^j + cxp ^-^fc p (d - g)j sin ~ 0) 
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y J n a 2 Vp; 4F 2 + 4^7-2 %o a 4a<5 (a 2 + <5 2 ) 2 V ' 



x < a 



(a 2 + 5 2 ) 2 + 2l3 2 7 2 (35 2 -a 2 )FF . , . 

^ '- 1 1 _ 2i/3J F - 2 7 2 F 

or + o z 



cxp ^-fc p ^ cos (—kp^j — exp ^— — k p (d — cos (^—k p (d — £)j 

OL ~\- 

cxp sin (~ fc P^) + CX P _ sin (^p(^ ~ ^) } ' 

As it follows from the expression (56), in regions I and II the longitudinal field 
ahead of bunch (£ > d) is missing. The transverse fields ahead of the bunch are 
independent of £ (i.e., are not modulated), are circularly polarized and proportional 
to the function F and thickness of the bunch d. Note that F is proportional to the 
difference between the group velocity of induced transverse fields (v g — k c 2 /uj ) 
and the bunch velocity. If the group velocity of the transverse wave coincides with 
the bunch velocity (F = 0), then the transverse fields ahead of the bunch are absent. 
Besides, the magnetic field ahead of the bunch is proportional to F and is absent 
in the long- wave (quasi-stationary) limit (fco = or F = 0). 

The transverse fields ahead of the bunch arise on account of the fact that the 
CPEM wave has a phase velocity (v v — w /fc > c ) exceeding the bunch velocity 
at any values of parameters of plasma and of the CPEM wave. Hence, some part 
of perturbations caused in plasma by the CPEM wave and bunch may have phase 
velocity higher than the bunch velocity and be even ahead of it. 

Inside and behind of bunch for values of parameters in regions I, II, III, as well 
as ahead of it in the region III (see the Eqs. (59), (61), (63), (64), (66), (67), (69), 
(70), (72)-(77)) the transverse fields are modulated and are circularly polarized. 
Though the polarization vector of the transverse wave describes a circle, its radius 
depends on the distance to the bunch £. Indeed, as it follows from Eqs. (29)-(31) 

E 12 + E 12 = E 2 (0 + E 2 {£) = EZuJS), (78) 

B' 2 + B' 2 = B 2 r {0 + B 2 {i) = B 2 max (0- (79) 

The amplitudes of electric (E max (£)) and magnetic (.B max (£)) fields are generally 
dependent on £. The modulated transverse waves in plasma arise in consequence 
of the excitation of two type waves, having frequencies ojlX± and wave vectors 
(lul/u)\±. Owing to the interaction of the pump wave with these induced waves, 
there arise waves with combined frequencies ujq — lul\±, luq + lul\± and combined 
wave vectors fc — (u)l/u)\±, fc + (u>l/u)X±, the interference between which gives 
a modulated wave. 

Here E r (£) or serve as a carrier wave if loo > <^lA± (in the dimcnsionless 

form a A > A±) and fc > (u>l/u)X± (in the dimensionless form /?\/a 2 A 2 — 1 > A±). 
Otherwise in expressions (29) and (30) the part of the carrier wave will play the 
functions cos£ and sin£. 

Below we shall analyze the obtained expressions in practically important case 
of the following values of parameters: no < 10 17 cm~ 3 , (uj p < 2 x 10 13 sec _1 ), Wl — 
10 18 W/cm 2 + 10 20 W/cm 2 (W L = cE$/4tt is the intensity of CPEM wave), lu ^ 
10 15 sec -1 , 7 ~ 10-hl0 3 . At the variation of CPEM wave intensity in the range from 



15 



10 18 W/cm 2 to 10 20 W/cm 2 the parameter a is changed in the interval a ~ 1.02-^2. 
For the parameter A we have: A > 50. 

In the region I 7 < A or 7 > A (the Eq. (40)). It is comparatively easy to 
analyze the Eqs. (56)-(64) in two particular cases of aA >7> 1 and 1 <C aA <C 7 
respectively. In the first case we obtain from Eqs. (34) and (35): 

4 7 4 F 2 * a 2 A 2 - » 7 2 » 1, (80) 

A^-l, A 2 h ~ 4 7 4 F 2 ~a 2 A 2 »A 2 _. (81) 

The transverse electric and magnetic fields ahead of the bunch are on the order 
of magnitude \£ + \ ~ \B + \ ~ E (rib/no)(k p d)j 2 /2 when a > 1. Inside and behind 
the bunch the amplitude of induced waves with frequency ujlX + is much less than 
the one for waves with frequency wz,A_. By the order of magnitude the amplitude 
of oscillations with frequency lj^X- in the longitudinal waves (Eqs. (59) and (62) 
is Eo(rib/no)a, i.e., is r y 2 (k p d)/2a times as less as the amplitude of transverse waves 
ahead of the bunch. For ^ £ ^ d and £ < the highest contributions to expressions 
(60), (61), (63), (64) are made by the first and third terms, the amplitude of the 
third term being approximately equal to E (nb/no)-f 2 a (i.e., j 2 times as much as 
the amplitude of longitudinal waves). Hence, in the region I S> E r (£) ~ 0, 

B r {0 > B i(0 - and E x S -£;(£) sin C, ^ = £ 4 (£)cosC, £ x S B r (£)cosC, 
B y = i?i(£)sin£. Besides, the oscillating terms in Eqs. (59)-(64) describe the 
waves, the length of which increases with the intensity of CPEM wave. Third, 
when the condition ird < a\ p is met (where A p = 2ix /k p is the wavelength of excited 
longitudinal waves in the absence of the CPEM wave), the oscillating terms in 
expressions (60), (61), (63), (64) exceed the first terms for narrow bunches. In case 
of wide bunches (ird > aX p ) the first term in Eqs. (63) and (64) exceeds the third 
term. Inside the wide bunch the analogous conclusion is valid close to the bunch 
boundaries when tt \d/2 — £| > aX p /2. Inside the wide bunch, in the vicinity of its 
center, tt \d/2 — £| < aX p /2, the main contribution to Eqs. (60) and (61) is made 
by the third terms. 

Behind the bunch, the amplitude of longitudinal and transverse waves is pro- 
portional to 2sin(7rd/aA p ) and is maximum at d — (n — l/2)aX p , where n = 1, 2, .... 
When the condition d = (aX p )n is fulfilled, the waves behind the bunch arc not 
excited. 

In Figs. 4-10 are shown the results of numerical calculations made for induced 
fields using Eqs. (56)-(64). The figures illustrate all the features discussed earlier. 
Thus, in the presence of a CPEM wave in the plasma, the one-dimensional bunch 
excites the waves with the wavelengths increasing as the intensity of CPEM wave. 
The dependence of longitudinal wave amplitude on the energy of bunch (on the 
rclativistic factor) is weak and grows with the intensity of CPEM wave, whereas 
the amplitude of the transverse wave is j 2 times as large as that of longitudinal 
wave. One can assert that at 7 3> 1 the induced wave is almost a transverse one. 

Now consider another limiting case with 7 3> aA »1. Here instead of Eqs. (80) 
and (81) we find 

4 7 4 ^7 2 (^) 2 »1, (82) 

A 2 _cl, A^7 2 (^) 2 »A 2 _. (83) 

When a > 1, the transverse fields ahead of the bunch are by the order of mag- 
nitude: \£ + \ ~ \B + \ ~ E {nb/n )(k p d)a 2 A 2 /2. Inside and behind the bunch the 
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amplitude of waves with frequency is again much less than the amplitude of 

waves with frequency wj,A_. The amplitude of oscillations with frequency in 
the longitudinal waves (Eqs. (59) and (62)) docs not change and is nearly equal to 
Eo(nb/no)a. In this case the magnitude of the longitudinal wave is by (k p d)aA 2 /2 
times less than the transverse wave amplitude ahead of the bunch. Inside and be- 
hind the bunch the primary role in the transverse waves is played again by the first 
and third terms. The amplitude of the third terms is E (nb/no)a(aA) 2 (i.e., by 
(aA) 2 times larger than the amplitude of longitudinal waves). Note that in the 
limit of very large values of 7 (7 » aA) the transverse fields are independent of 
it, but the dependence on intensity of CPEM wave becomes more salient. All the 
remaining features mentioned above for the case of 7 <C aA are still valid here. 

As we saw in Sec. Ill, 7 ~ aA in the second region of values of a, A, 7. In 
this case the function F in Eqs. (56)-(58), (65)-(70) takes on the values F± (see 
the expression (38)), the region II being determined by expressions (36) for 7 3> 1 
and an arbitrary value of parameter a, and by expressions (37) for 7 3> 1 and 
1 < a < 00(7) (where 00(7) is determined by the expression (39)). The values 
F = ±F + correspond to Eq. (36), and the values F = ±F_ correspond to Eq. (37). 
Respectively, in the expression (45) one is to take the sign "+" if F = ±F + , and 
the sign "— " if F = ±F_. From expression (45) for values 7 > 1, A > 1 we have 
A ~ 7 1/2 > 1 if a > 1 + l/4 7 2 , and A < 1 if 1 < a < 00(7) ^ 1 + I/I67 2 . In 
the first case the wavelength of the waves excited inside and behind the bunch is 
decreased as 7 increases and is much less than the wavelength in the second case, 
when that is on the order of aA p . 

As it follows from Eqs. (65)-(70), at the distance from the bunch the amplitudes 
of waves to be excited increase proportional to the distance to bunch and beginning 
from some value £ = £0 the fields may be of the same order of magnitude as the 
amplitude of CPEM wave. Surely, in case of £ > £0 the Eqs. (68)-(70) are not 
adequate for the description of induced fields and for their treatment one is to avail 
of the nonlinear approach. 

The obtained results for region II were further analyzed numerically. In Figs. 
11-14 the coordinate dependence of induced longitudinal and transverse fields for 
the following values of CPEM wave, plasma and bunch parameters is shown: n& = 
10 14 cm" 3 (Fig. 11), n b = 10 13 cur 3 (Figs. 12-14), n a = 10 17 cur 3 , lo q = 3.77 x 
10 15 sec-\ 7 = 210, d = 2k' 1 , k p 1 = 1.7 x l(T 3 cm. The strength of CPEM wave 
field is varied from the value E = 8.5 x 10 8 V/cm to E Q = 1.85 x 10 10 V/cm. The 
^-dependence of i? max in Figs. 11-14 is not shown, since the difference of B max from 
-Emax(C) is insignificant as it follows from Eqs. (66), (67) and (69), (70) for aA » 1, 
7 3> 1. From Fig. 12 one can estimate the distance £0, at which the amplitude of 
transverse waves is of the order of CPEM wave amplitude and the linear approach 
is of no sense any more. From Fig. 12 one gets the following estimate: £0 — 60fc~ x . 
So, the linear treatment is valid at distances of several wavelengths (2ir/k p ). 

In the region III the fields are described by Eqs. (71)-(77). At large distances 
from the bunch |£| > a(<5fc p ) -1 , |£ — d\ > a(£fc p ) _1 , the longitudinal waves are 
attenuated, and transverse waves are described by the first terms in Eqs. (72), (73) 
and (76), (77) that are independent of £ (i.e., are not modulated). The ^-dependence 
of transverse and longitudinal waves is given in Figs. 15-21. It follows from these 
figures that the longitudinal field and amplitudes E r , Ei, B r are antisymmetric with 
respect to the bunch center (£ = d/2 plane), and the amplitude Bi is symmetric 
with respect to it. Besides, at large distances the phase of transverse waves ahead 
of the bunch differs from the phase of transverse wave behind the bunch by n, and 
the transverse oscillations behind and ahead of the bunch are in the counterphasc. 
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5 Summary 



The purpose of this work was to investigate the excitation of EM wake waves in 
electron plasma by one-dimensional bunch of charged particles in the presence of in- 
tense CPEM wave. The interaction of the pump wave with a plasma was described 
by means of Maxwell equations and relativistic nonlinear hydrodynamic equations 
of plasma. Then, we have considered small perturbations of plasma due to the 
presence of one-dimensional bunch. A general expressions obtained for EM wake 
field components was analyzed in three particular ranges of parameters of the pump 
wave, bunch and plasma. In the first range of parameters the amplitude of trans- 
verse components of induced waves is shown to grow as the bunch energy and after 
some definite value of the relativistic factor of the bunch to be almost independent 
of the energy and increase proportional to the intensity and frequency of the CPEM 
wave. The dependence of the longitudinal component of induced waves on the rel- 
ativistic factor of the bunch is weak. Its amplitude and wavelength grow as the 
intensity of CPEM wave. In the second range of parameters the amplitude of the 
wave excited by the bunch is a linear function of the distance to the bunch. In the 
third range of parameters the longitudinal component of induced fields is localized 
near the bunch boundaries and exponentially decreases with increasing distance 
from these boundaries. The amplitudes of transverse components of induced waves 
reach a constant value with the distance from the bunch boundaries. In general, 
the transverse waves in all three mentioned ranges of values of the CPEM wave, 
plasma and bunch parameters are modulated and are circularly polarized with the 
same type (right-hand or left-hand) of polarization as the pump wave. 

In conclusion it is worthwhile to make two remarks. First, although at the 
consideration of EM wake wave generation we took the external CPEM wave to 
be strong, it was assumed monochromatic. The intensity of present-day sources of 
monochromatic EM waves does not exceed 10 16 W/cm 2 and the values of amplitudes 
used for numerical estimates and calculations arc much higher than this value. In the 
nearest future we plan to study the excitation of EM wake waves by one-dimensional 
bunch in the presence of high power laser pulses with the intensity ~10 20 W/cm 2 . 
The results of this study will be presented in a forthcoming work. Second, it follows 
from Eqs. (10)-(14) that in the plasma there arise parametrically coupled linear 
waves. In case of high bunch densities (rib > no) the linear treatment is of nosense 
any more and it is necessary to examine explicit system of nonlinear Eqs. (l)-(5). 
This issue will also be studied in future work. 
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6 Appendix 



In Section II a set of inhomogeneous differential equations with constant coefficients 
for induced variables £ ± , , V ± , E z and n (see Eqs. (16)-(20)) has been found. 
The solution of this set of equations was obtained with the help of Green's function. 
Below we shall briefly discuss another method for solution of the set of Eqs. (16)- 
(20), i.e., we shall obtain a differential equation for one of the induced quantities, 
i.e., for the density of induced charge n . The equations for other variables will be 
obtained in an analogous way. 

Now consider the Eq. (20) and find the differential operator the action of which 
on B + + B~ and V + + V~ give a function of n . Let us introduce the operators 

After successive action of operators Q~ and Q + on both the sides of the Eq. (18) 
we obtain the equation that comprises only the functions V ± (^) and n (£): 



i -f) • (v «)±. 



n (0 \ , A : 



n 



«(0 



n 
(A-2) 

Now,Jntroduce another operator S, that is symmetric with respect to operators Q~ 
and Q + , and indeed: 



S = Q + Q- £ (l - f ) (Q- + Q + ) + ^ (1 - Pi) . (A. 



3) 



After some transformations one can write the Eq. (A. 2) with the help of operator 
S 



SV ± (Z)=±iv, 



■ n 1 



1 



cr 



(1 ~ PD 



no 



(A.4) 



From Eqs. (A.4) and (17) the required equation follows for variables V + V + and 
Br + B+: 



S[V"(0+V + (0]=-4«e^(fco-^) 



wo\ g n (g) 
c / dl; no 



S [B-(0 + B+(0] = -^eno(3 e j (q~ + 



9^ 



+ ifc Q + - Q~ 



(A.5) 
5 I n (0 



c 2 <9£ J no 
(A.6) 

Now, acting by the operator S on both the sides of the Eq. (20) and using the rela- 
tions (A.5), (A.6), we finally arrive at the equation that contains only the variable 

n 



d 4 * d2 o 

w +A w +B 



no 



u 2 \oe 



+ c 



n 



(A.7) 



where 
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A=^(l-Pl)+^(k,-^)\ (A.8) 
c z V c 



^ = 7 2 



c 



(A.9) 
(A.10) 



One can obtain the equations for variables £ ± , B ± , V ± , E z in an analogous way. 

Prior to the consideration of Eq. (A. 7) note first of all that in the absence of 
CPEM wave {fi e = 0) one has the following relation from Eqs. (A.7)-(A.10): 

* +A * +B ).(* +c )(* + 4) «A.n, 



and Eq. (A. 7) passes into the well known equation for the induced charge obtained 
in Rcfs. [24, 25] 

( d 2 L0 2 \ , LO 2 

{w + *) n ® = -* Nb ®- (A ' 12) 

In the limit of strong CPEM wave ((3 e ~ 1) we have from expressions (A.7)-(A.10) 

A~C ~4%1 + Pr 2 , B~0, (A.13) 
& 

8 2 ( 8 2 \ , 

^(^2+<?o)n(0 = 0, (A.14) 

where qo = \J~A = 2(wo/c)(l + /3) _1 . It follows from expression (A.14) that in this 
case the equation for n will be a homogeneous differential equation. Analogous 
equations are obtained for the rest of the induced variables. So, in the presence 
of high intensity CPEM wave the bunch will not altogether perturb the state of 
homogeneous plasma established by the external wave and all induced variables 
will be zeros. 

For an external wave of arbitrary intensity a non-homogeneous equation of the 
fourth degree is obtained for n (£) (Eq. (A. 7)), the coefficients of which depend 
on the intensity of CPEM wave. The characteristic equation corresponding to Eq. 
(A. 7) establishes in general case a dispersion law for two types of induced waves 
that are parametrically coupled due to the presence of the external wave (it is 
easy to see that it coincides with Eq. (33)). In case of A 2 > 4B Eq. (A. 7) 
describes the oscillations of induced charge density with frequencies of lul\- and 
ujl\ + respectively. This type of solution corresponds to the first range of parameter 
values that was studied above. In case of A 2 = AB the characteristic equation has 
multiple solutions and respectively the solutions of Eq. (A. 7) describe the waves, 
the amplitude of which increases as £ (range II). At A 2 < AB the solutions of the 
characteristic equation are complex and the waves excited by the electron bunch 
exponentially decrease with the distance from the bunch (region III). 

To obtain a single- valued solution of the problem one is to supplement Eq. (A. 7) 
with boundary conditions. As such in case of an arbitrary, smoothly changing 
density profile of the bunch -/Vb(£) when £ — ► ±oo, one can use the equality of 
function n (£) and of its first three derivatives at £ — > +oo (or £ — ► — oo) to zero. 
A distinguishing feature of the one-dimensional bunch with sharp boundaries (the 
profile of this kind for iV&(£) we have used in calculations of induced EM fields) is 
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the fact that the even derivatives of the function n (£) are jump functions at the 
boundaries of the bunch. Indeed, n (£) is continuous at the boundaries of the bunch 
(for £ = and £ = d). As it follows from Eq. (20), dn /d£ is also continuous on the 
boundaries of the bunch. We shall obtain by integrating Eq. (A. 7) over the variable 
£ in an infinitesimal range in the vicinity of points £ = or £ = d, that d 3 n /<9£ 3 
is also continuous on the boundaries of the bunch. Now consider Eq. (A. 7) in the 
points £ = £o — and £ = £ + (where £ = or £ = d) respectively. The 
subtraction of the obtained relations gives the following condition: 

Co+o 2 

= -^Can b . (A.15) 

Co-o 

where a — +1 for £o = and a = — 1 for £o = d- Besides the condition (A.15) 
and the conditions of continuity of functions of n , dn /d£ and d 3 n /<9£ 3 at the 
boundaries of the bunch, there is another boundary condition on the induced vari- 
ables describing the longitudinal waves in plasma (i.e., on n (£), v z {£) and E z (£)). 
According to this condition, owing to the equality of phase velocities of longitudinal 
waves induced in the plasma to the velocity of bunch ahead of these, the functions 
n (£) (and, therefore, v z (£)) and E z {£) turn zero. Thus, all the mentioned boundary 
conditions for a bunch with sharp boundaries together with Eq. (A. 7) specify the 
unique solution of Eq. (A. 7). The equations and boundary conditions for transverse 
induced variables {S^, B ± , V ± ) the solutions for which coincide with expressions 
(22)-(25), may be obtained in an analogous way. 




References 

[1] E. Esarcy, P. Spranglc, J. Krall, and A. Ting, IEEE Trans. Plasma Sci. 24, 
252 (1996). 

[2] Ya. B. Fainberg, Fiz. Plazmy 23, 275 (1997) [Plasma Phys. Rep. 23, 251 
(1997)]. 

[3] A. Ts. Amatuni, S. S. Elbakian, A. G. Khachatryan, and E.V. Sekhpossian, 
Part. Accel. 51, 1 (1995). 

[4] P. Chen, Part. Accel. 20, 171 (1987). 

[5] P. Chen, J. M. Dawson, R. W. Huff, and T. Katsouleas, Phys. Rev. Lett. 54, 
693 (1985). 

[6] R. D. Ruth, A. W. Chao, P. L. Morton, and P. B. Wilson, Part. Accel. 17, 
171 (1985). 

[7] R. Keinigs and M. E. Jones, Phys. Fluids 30, 252 (1987). 

[8] A. Ts. Amatuni, E. V. Sekhpossian, A. G. Khachatryan, and S. S. Elbakian, 
Fiz. Plazmy 21, 1000 (1995) [Sov. J. Plasma Phys. 21, 945 (1995)]. 

[9] A. G. Khachatryan, A. Ts. Amatuni, E. V. Sekhpossian, and S. S. Elbakian, 
Fiz. Plazmy 22, 638 (1996) [Sov. J. Plasma Phys. 22, 576 (1996)]. 



21 



[10] Ya. B. Fainberg, N. Ayzatskij, V. Balakirevet al., Proc. of XV Int. Workshop 
on Charged Part. Linear Accel., 1997, Alushta, p. 16 (in Russian). 

[11] S. S. Elbakian, E. V. Sckhpossian, and A. G. Khachatryan, Preprint YerPhl- 
1511(ll)-98 (E-Print: ps/9804017). 

[12] A. Ts. Amatuni, E. V. Sekhpossian, and S. S. Elbakian, Fiz. Plazmy 12, 
1145 (1986). 

[13] J. B. Rosenzweig, Phys. Rev. Lett. 58, 555 (1987). 

[14] J. B. Rosenzweig, IEEE Trans. Plasma Sci. 15, 186 (1987). 

[15] J. B. Rosenzweig, Phys. Rev. A 38, 3634 (1988). 

[16] A. G. Khachatryan, Phys. Plasmas 4, 4136 (1997). 

[17] J. M. Dawson, Phys. Rev. 133, 383 (1959). 

[18] A. I. Akhiezer and R. V. Polovin, Zh. Eksp. Teor. Fiz. 30, 915 (1956) [Sov. 
Phys. JETP 3, 696 (1956)]. 

[19] N. E. Andrccv, L. M. Gorbunov, and R. R. Ramazashvili, Fiz. Plazmy 23, 
303 (1997) [Plasma Phys. Rep. 23, 277 (1997)]. 

[20] A. M. Kalmykov and N. Ya. Kotsarcnko, Izvestya VUZ 19, 1481 (1976) (in 
Russian) . 

[21] L. Stenflo, Plasma Physics 19, 1187 (1977). 

[22] C. Max and F. Perkins, Phys. Rev. Lett. 29, 1731 (1972). 

[23] C. S. Liu and V. K. Tripathi, Interaction of Electromagnetic Waves with Elec- 
tron Beams and Plasmas (World Scientific Publishing, Singapore, 1994). 

[24] M. E. Jones and R. Kcinigs, IEEE Trans. Plasma Sci. 15, 203 (1987). 

[25] Hyun-Soo Kim, S. Yi, A. Amin, and K. E. Lonngren, Phys. Rev. E 50, 3962 
(1994). 



22 



FIGURE CAPTIONS 



FIG. 1. The line A = A (a), where the expression under the square root in Eq. 
(34) turns zero. The region III is contained within the lines, the region I is outside 
the lines, and the curve A = A(a) is the region II. The plot is made for 7 = 1.5. 

FIG. 2. The line A = A(a) for 7 = 20. 

FIG. 3. The line A = A(a) for 7 = 100. 

FIG. 4. The ^-dependence of the induced longitudinal electric field in the region 

I of the values of parameters. The field is measured in the units of 10 5 V/cm, £ is 
measured in units of k~ x = 1.7 x 10~ 3 cm (zq = k p ^). The curves were calculated for 
the following values of parameters: no = 10 17 cm~ 3 , nb = 10 14 cm~ 3 , uoq = 3.77 x 
10 15 sec -1 , 7 = 50, k p d = 20. The dotted line corresponds to the absence of CPEM 
wave (Eq = 0), the dashed line corresponds to the value Eq = 1.3 x 10 1:L V/cm, the 
solid line corresponds to the value E — 2.5 x 10 11 V/cm. 

FIG. 5. The ^-dependence of maximum value of induced transverse electric field 
^£ mra = a/-E 2 (£) + i? 2 (£)^in the region I. i? max is measured in units of 10 9 V/cm, 

£ in units of 1.7 x 10 -3 cm. The thickness of the bunch is lOfc" 1 . The dotted line 
corresponds to E = 6.7 x 10 10 V/cm, the dashed line to Eq = 8.8 x 10 10 V/cm, 
the solid line to Eq = 1.09 x 10 11 V/cm. The remaining parameters coincide with 
parameters given in Fig. 4. 

FIG. 6. The ^-dependence of field E r (£) (in units of 10 5 V/cm, and £ in units of 
kp 1 = 1.7 x 10~ 3 cm). The parameters and notations are the same as in Fig. 5. 

FIG. 7. The ^-dependence of field Ei(£) (in units of 10 9 V/cm, and £ in units of 
kp 1 = 1.7 x 10~ 3 cm). The parameters and notations are the same as in Fig. 5. 

FIG. 8. The ^-dependence of the maximum value of induced magnetic field 

(5 max = y/Bt(Z) (in units of 10 9 V/cm, and in units of k' 1 = 1.7 x 

10~ 3 cm). The parameters and notations are the same as in Fig. 5. 

FIG. 9. The ^-dependence of magnetic field ( m units of 10 9 V/cm, and £ 

in units of k~ x = 1.7 x 10~ 3 cm). The parameters and notations are the same as in 
Fig. 5. 

FIG. 10. The ^-dependence of (in units of 10 5 V/cm, and £ in units of 

kp 1 = 1.7 x 10~ 3 cm). The parameters and notations are the same as in Fig. 5. 
FIG. 11. The ^-dependence of induced longitudinal electric field in the region 

II (in units of 10 6 V/cm, and £ in units of k~ 1 = 1.7 x 10~ 3 cm). The curves 
were obtained for values of parameters: n b — 10 14 cm~ 3 , n = 10 17 cm~ 3 , ui a = 
3.77 x 10 15 sec _1 , 7 = 210, k p d = 2. The dotted line corresponds to the value 
Eq = 8.5 x 10 8 V/cm, the solid line corresponds to the value Eq = 1.85 x 10 10 V/cm. 

FIG. 12. The ^-dependence of E max in the region II (in units of 10 9 V/cm, and 
£ in units of A;" 1 = 1.7 x 10~ 3 cm) for the values of parameters: rib — 10 13 cm~ 3 , 
Eq = 1.85 x 10 10 V/cm. The remaining parameters are the same as in Fig. 11. 

FIG. 13. The ^-dependence of fields £y(£), B r (0 in the region II (in units of 
10 9 V/cm, and £ in units of k^ 1 = 1.7 x 10~ 3 cm) for the values of parameters: 
rib = 10 13 cm~ 3 , kpd = 2. The values of remaining parameters are the same as in 
Fig. 11. The dotted line corresponds to the value of E r (£), the solid line corresponds 
to B r (0- 

FIG. 14. The ^-dependence of fields Sj(£), in the region II (in units of 

10 9 V/cm, and £ in units of fc" 1 = 1.7 x 10~ 3 cm). The parameters are the same 
as in Fig. 13. The dotted line corresponds to Ei(£), the solid line corresponds to 

FIG. 15. The ^-dependence of the induced longitudinal field in the region III (in 
units of 10 2 V/cm, and £ in units of fc^ 1 = 1.7 x 10~ 3 cm). The parameters are: rib — 
10 14 cm- 3 , riQ = 10 17 cnT 3 , cu = 3.77 x 10 15 sec _1 , 7 = 210, k p d = 5. The dotted 
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line corresponds to the value E = 4.09 x 10 9 V/cm, the dashed line corresponds to 
E Q = 1.16 x 10 10 V/cm, the solid line corresponds to E Q = 1.6 x 10 10 V/cm. 

FIG. 16. The ^-dependence of E max in the region III (in units of 10 8 V/cm, and 
£ in units of fc" 1 = 1.7 x 10~ 3 cm) for rib = 10 13 cm~ 3 . The remaining parameters 
and notations are the same as given in Fig. 15. 

FIG. 17. The ^-dependence of the field E r (£) in region III (in units of 10 8 V/cm, 
and £ in units of k~ x = 1.7x 10~ 3 cm) for n& = 10 13 cm~ 3 . The remaining parameters 
and notations are the same as given in Fig. 15. 

FIG. 18. The ^-dependence of the field Ei(£) in region III (in units of 10 8 V/cm, 
and in units of k^ 1 = 1.7 x 10 _3 cm) for rib = 10 13 cm -3 . The remaining parameters 
and notations are the same as given in Fig. 15. 

FIG. 19. The field -B max (£) in the region III. The values of parameters, mea- 
surement units and notations are the same as in Fig. 16. 

FIG. 20. The field B r (£) in the region III. The values of parameters, measure- 
ment units and notations are the same as in Fig. 16. 

FIG. 21. The field £?i(£) in the region III. The values of parameters, measure- 
ment units and notations are the same as in Fig. 16. 
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